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Abstract. 

The energy spectrum of cake shape normal - superconducting systems is calculated 
by solving the Bogoliubov-de Gennes equation. We take into account the mismatch 
in the effective masses and Fermi energies of the normal and superconducting regions 
as well as the potential barrier at the interface. In the case of a perfect interface 
and without mismatch, the energy levels are treated by semi-classics. Analytical 
expressions for the density of states and its integral, the step function, are derived and 
compared with that obtained from exact numerics. We find a very good agreement 
between the two calculations. It is shown that the spectrum possesses an energy gap 
and the density of states is singular at the edge of the gap. The effect of the mismatch 
and the potential barrier on the gap is also investigated. 
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1. Introduction 

At the interface of the mesoscopic hybrid normal - superconducting (NS) junctions the 
underlying physics is controlled by the coherent evolution of electrons into holes due 
to the mechanism called Andreev reflection pQ. The growing interest in studying the 
role of the phase coherent Andreev reflection is motivated by the recent technological 
progress in manufacturing almost ballistic semiconductors of mesoscopic size coupled to 
a superconductor (for overview of the recent progress in this field see e.g. Ref. [21III1IUE])- 
A ballistic normal dot weakly coupled to a superconductor is commonly called an 
Andreev billiard. The excitation spectrum (Andreev states) of such NS systems can 
be calculated from the Bogoliubov-de Gennes equation t 6j (BdG). A film of normal 
metal in contact with semi-infinite superconductor was considered by P. G. de Gennes 
and D. Saint- James can be regarded as the first such Andreev billiard. The 
bound states of Andreev billiards have been extensively studied over the past ten 
years |Hl UHI [TTJ [T21 [THl U^l USl HSl UHl [THl ED! - The excitation spectrum were studied 
for SNS junctions [21112211231211123, too. The energy levels of the strong magnetic field 
were investigated for a semi-infinite N region in contact with a semi-infinite S region j2B] 
and ring-shaped Andreev billiards [27 j . The Andreev bound states for superconducting- 
ferromagnetic systems have also been studied (see e.g., Ref. [2E] and refs. therein). 

An important physical quantity for Andreev billiards is the density of states (DOS). 
For perfect NS interface (in the absence of tunnel barrier at the NS interface, and the 
mismatch in the Fermi energies and effective masses between the N and S regions) a 
semi-classical Bohr-Sommerfeld approximation of the DOS was derived for excitation 
energies close to the Fermi energy CHI El El El- R was shown that the DOS can 
be related to the purely geometry-dependent path length distribution P(s) which is the 
classical probability that an electron entering the N region at the NS interface returns to 
the interface after a path length s. In recent works PHI 1201 EHj an improved semi-classical 
Bohr-Sommerfeld approximation of the DOS has been presented, in which the energy 
dependent phase shift for Andreev reflections is properly taken into account. It has been 
shown that the exact quantum mechanical calculations for different Andreev billiards 
with perfect NS interface gives good agreement with the semi-classical predictions for 
all energies below the bulk superconducting gap (not only at energies close to the Fermi 
energy). For a NS system possessing a finite cutoff in P(s), a gap appears in the energy 
spectrum which can be as large as half of the bulk superconducting gap jTHj. Based on 
the semi-classical expression of the DOS, a simple formula for this gap was also derived 
for such Andreev billiards. 

In this paper our aim is twofold. First, we propose a new class of Andreev billiards 
which can also exhibit a gap owing to the finite cutoff in P(s). Second, we study how 
the gap is suppressed for realistic NS systems in which the NS interface is non-perfect. 
To this end we consider a 'cake shape Andreev billiard' shown in Fig. ^ in which the 
boundaries of the N and S regions are formed by two radii at an angle a and the NS 
interface is an arc between them. 
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Figure 1. A cake shape Andreev billiard formed from a cake shape of normal dot (N) 
in contact with a superconductor (S). 



Recently Adagideli et al. [16] investigated the Andreev billiards formed from an 
arbitrary shape of normal region surrounded by superconductors. Their theoretical 
predictions have been tested numerically in the case of circular shape of normal region. 
Similarly, the circular case has been studied by Stone |T7]. This NS system is very 
similar to our cake shape Andreev billiards for a = 2n. The only difference is that in 
our case there is an infinite potential barrier along a radius. 

The text is organized as follows. In Sec. 121 a quantization condition (secular 
equation) is derived from the BdG equation for the general case of non-perfect NS 
interface. In Sec. |3] the density of states is derived from semiclassics for the case of 
perfect NS interface. Our numerical results for the gaps in case of non-perfect NS 
interface are presented in Sees. |^ and [31 Finally, the conclusions are given in Sec. |3 

2. Exact quantum calculation 

In this section we derive a secular equation for the energy levels of cake shape Andreev 
billiards shown in Fig. ^ in the general case of non-perfect NS interface. The BdG 
equation for the NS systems can be written as 



where H = p 2 /2m e g + V(r) — /x is the single-particle Hamiltonian, fi = E F , E F are 
the Fermi energies, m e g = mN,ms are the effective masses in the N/S regions, \& is a 
two-component wave function, E is the quasi-particle energy measured from the Fermi 
energy E F N \ The tunnel barrier V(r) at the NS interface in polar coordinates (r,(p) is 
modelled in a usual way by V(r, ip) = Uq S(r — R). We also adopt the usual step-function 
model for the pair potential and take A(r) = AO(r — R). An infinite potential is 
assumed at the straight segments of the boundary of the NS system, i.e., the Dirichlet 
boundary conditions are applied for the wave function ^. In this work we are interested 
in the discrete energy levels, i.e., it is assumed that < E < A. 

It is easy to see that the Hamiltonian is separable in polar coordinates. The ip 
dependence of wave function \1/ is identical in the N and S regions. After separating the 




(1) 
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ip dependence of \I/ the BdG equation reduces to two Bessel equations [30] for the radial 
dependence of the two components of ^ . Thus, one can show that the ansatz for the 
wave functions in the N region can be written as 

a +Ju m (k + r 

while in the S region the wave function has the form 



sin (v m ip) 



(2) 



where v m = ^ (m = 1,2, 
functions [HO], and 









(r) 



sin (v m (p) , (3) 



H^{q + r) + c. 
J Vm {x) and H^' 2 \x) are the Bessel and the Hankel 
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(4) 
(5) 



3 ±i arccos(-B/A) 



7± = e ~" v ~'~'- ( 6 ) 
Here +/— refer to the electron/hole like quasi-particle excitation and the Fermi wave 
numbers in the N/S regions are given by k^ = \J 2m^E^ /Ti 2 and k^ = J 2ms-E' F S ^ /h 2 . 
The eigenfunctions of the energy levels are labeled by a fixed integer number m. In the 
N region only the Bessel function can be used for the radial wave function since the 
Neumann function is singular at the origin. The wave function \1/ in the S region must 
tend to zero as r — > oo. This condition can be satisfied by choosing the appropriate 
Hankel function for the electronic and hole like component of ^ in the following way. 
It is known [30 j that for large r the Hankel functions H^(q+r) ~ y / 2/(g + r) e iq+r . Since 
q + has a positive imaginary part one finds that Hff'(q + r) — ► for r — > oo. Similarly, 
H^(q_r) — > for r — > oo. 

The 99 dependent part of the wave functions in the N and S regions ensures that 
they satisfy the Dirichlet boundary conditions at the straight segments of the Andreev 
billiard. At the NS interface the matching conditions [3*T] ITH] should be applied. The 
four coefficients a±, c± in Eqs. (J2J) and (jUJ) are determined from the matching conditions: 



m 
d 

dr 



^(N) _ ^^(S) 
m m 



2m* 



WL S) 



-R 



(7) 
(8) 



r=R h 

Substituting the ansatz given by Eqs. (J2J) and Q into Eq. (J2J) one finds the following 
secular equation for the eigenvalues E of the NS system for fixed mode index m: 

Im{ 7+ D«(E)D«(E)} = 0, (9) 
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where Im{.} stands for the imaginary part and the two by two determinants Dffl(E) 
are given by 



r>(+) r>(+) 

-^11 ^12 

D (+) n(+) 

-^21 -^22 



(10) 



D { m \E)=[D^\-E)]\ (11) 

with matrix elements 

D[V = J„ m (k + R), (12) 

D& = HM(q + R), (13) 

M^^UM?), (14) 

* =" 2k^ZHl]}( q+ R) + ^ q+ H^ {q+ R), (15) 

m s 

and Z = Uokp 1 ^ /(2E^) is the normalized barrier strength used by Blonder et al. |29j . 
The primes denote the derivatives of the Bessel and Hankel functions with respect to 
their arguments. For a given quantum number m the energy levels of cake shape Andreev 
billiards for non-perfect interface are the solutions of the secular equation given by 
Eq. Q. This equation is exact in the sense that the usual Andreev approximation 
is not used. In the Andreev approximation it is assumed that A / E F <C 1 and 
quasiparticles whose incident /reflected directions are approximately perpendicular to 
the NS interface [3]. 

3. Perfect NS interface 

In the case of a perfect NS interface, the semi-classical expression for the DOS can be 
derived from the secular equation (|9|) . We assume that there is no mismatch and tunnel 
barrier at the NS interface (wn = ms, E F N ^ = E F ^ and Z = 0). For simplicity, we shall 
omit the superscript N and S in the wave numbers and the Fermi energies. 

In the Andreev approximation (A/E F <C 1) we take k + ~ q + in places where they 
are multiplied by the Bessel or Hankel functions in D^'(E) given by Eq. (jlU|) but in 
the arguments of the Bessel and Hankel functions we keep them to be different. To 
approximate D^\E) the Debye asymptotic expansion [^U] of the Bessel functions will 
be used for \v m \ < k + R — ^/k + R. For Hankel functions of a complex argument z the 



principal asymptotic forms jHOl HjP(z) ~ y2/(nz)e l ^ z 2U7T *> for \z\ — > oo will be 
applied. After some algebra, the quantization condition (jUJ) can be simplified as 



$(k + R, u m ) — $(k_R, v m ) — arccos 4 



7T 

where 



n, (16) 



$(x, u) = y/x 2 — v 2 — \u\ arccos — , (17) 

x 

and n is an non-negative integer. This is the Bohr-Sommerfeld quantization condition 
for cake shape Andreev billiards. The function d(k±R 1 v m ) is proportional to the radial 
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action of the electron/hole like quasiparticles \32\ . For a given m and n the solution of 
(fTSJl gives the E mn energy level in Bohr-Sommerfeld approximation. 

We now introduce the so-called step function N(E) which is the number of energy 
levels below the energy E. This is indeed the integrated DOS g{E) of the Andreev 
billiards: 

N(E) = [ E q(E') dE' = £ 8(£ - E mn), (18) 

where Q(x) is the Heaviside function. To calculate the step function in Bohr-Sommerfeld 
approximation, Eq. (jl6j) is Taylor expanded in terms of the small quantity E/E-p. 
Keeping the leading terms one finds 

E I 2 E 

— y (kpR) — = nn + arccos — . (19) 
E-p "' A 

Substituting the solutions E mn of this equation into (fTSj) gives the step function Nbs(E) 

in Bohr-Sommerfeld approximation. Applying the Poisson summation formula [SSI E2] 

in the sum over m of N(E) and keeping only the non-oscillating term we find 

N BS (E) = £ / a &(E- E mn ) dm=J2 m*(E, n), (20) 

n=0 2 n=0 

where m*(E, n) is the solution of Eq. (|T9*j) for m at a given E and n, and M = - kpR is 
the largest m for which the argument of the square root in the left hand side of (|19|) is 
still positive. 

After a simple algebra we obtain the final form of the step function in Bohr- 
Sommerfeld approximation: 



N BS (E) = MY, {l-F[s n (E)]}, (21) 



ra=0 



where 



F(s) 



I* P(s') ds', (22) 
J o 



P(s) =-—2 , 2 &{2R-a), (23) 



, , n7r + arccos x / x 

«»(^) = E/A ^ (24) 

M = - kpR and ^ c = hv-p/(irA) = 2Ep / '(irkpA) is the coherence length. It can be shown 
that P(s) is the classical probability that an electron entering the billiard at the NS 
interface returns to the interface after a path length s. One can see that P(s) depends 
only on the radius R of the N region but it is independent of the angle a. This feature 
can be understood from a classical point of view. The unfolded trajectory of the particle 
starting and ending at the NS interface shown in Fig. El is always a chord of a circle of 
radius R, thus its length is independent of the angle a. The path length distribution 
P{s) is normalized to one, i.e., P{s) ds = 1. The integrated path length distribution 



is F{s) = [1-Jl- (s/2RY Q(2R - s). 
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Figure 2. Unfolded trajectory in the cake shape Andreev billiard. 



The DOS Qbs(E) = dN-Bs(E)/dE in Bohr-Sommerfeld approximation can easily 
obtained from (j21j) and reads as 

s 1 



roc 

g BS (E) = M / dsP(s 
Jo 



X 

71=0 



hv F VA 2 - E 2 

«j f s E ( E\\ 

5 ( \rm + arccos — ] ) , (25) 

^ \hv F V A// v ; 

where v p is the Fermi velocity. Note that expressions (|21|) and (|23j) for the step function 
and the DOS in the semi-classical approximation are the same as that presented in 
Ref. [TTJ] ED] for general shapes of Andreev billiards. Our formula for the DOS, valid 
for all energy E < A, can be regarded as an extension of that derived by Melsen et 
al. [01110], Lodder and Nazarov [TH], Schomerus and Beenakker [13], and Ihra et al. [TH] 
in the limit £ < A. In these works, the energy dependent phase shift — arccos(i?/A), 
due to Andreev reflections, was approximated by 7r/2, while in our result this phase 
shift is fully incorporated. 

In Fig. [S] the numerically calculated exact step function N(E) obtained by solving 
Eq. (J0J) together with its semi-classical approximation Nbs(E) from Eq. (|2*T|) as functions 
of the energy E below the bulk gap A are shown for angles a = ir, tt/2,tt/3. One can 
see that the agreement between the two results is excellent for all energy E < A. To 
see the difference between the two calculations, the inset shows an enlarged part of the 
main frame. 

The step functions shown in Fig. E] have cusps at some energies. This implies that 
the DOS has peaks at these energy values. Since the semi-classical approximation of the 
exact step function is very good one can derive an expression for the positions of these 
peaks starting from (j21j) . It is clear from (j2*H|) that P(s) is singular at s — 2R. This 
implies that the DOS will be singular at those energies when s n (E) = 2R. Using (|24j) 
and the approximation arccos (E/A) m tt/2 — E/A the positions of these singularites 
in the DOS is given by 

4 Sins) (n + l/2)7T 
A 1 + 2R/(^ C ) {0) 

valid for all integral n for which E^ mg ^ < A. It is worth mentioning that the positions 
of the singularities are independent of the angle a of the cake shape Andreev billiards. 
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E/A 

Figure 3. The exact step function N(E) (solid lines) from Eq. © and its semi- 
classical approximation N-qs(E) (dashed lines) given by Eq. (|21|l for the cake shape 
Andreev billiards of angles a = ir, 7r/2,7r/3 as functions of E/A. The inset shows an 
enlarged part of the main frame. The parameters are k^R = 350 and A/Ep — 0.03. 
With these data £ C /-R = 0.06. The numerically obtained exact gap from Eq. is 
.Egap/A = 0.1355, while its semi-classical value from Eq. (|2Tj) is E gap /A = 0.1366. 



Figure El also shows that in the energy spectrum of the system there is a quite 
large gap. This gap always exists for cake shape Andreev billiards. To derive a formula 
for the value of the gap one can also start from the semi-classical approximation. It is 
obvious that P(s) possesses an upper cut-off, i.e., P(s) = if s > s max = 2R. Thus, from 
Eq. ()24|l it follows that the energy spectrum has a lower bound. An approximate value of 
this energy can be obtained by setting n = in ()26j) since the position of the singularity 
and the upper cut-off of P(s) is identical for cake shape Andreev billiards. The lowest 
energy level (the gap of the spectrum) in semi-classical approximation becomes 

A l + 2R/(^ c y 1 ] 

For £ c -C R the gap is E gap /A ps 7r 2 ^ c /(4_R). It is clear that the value of the gap is also 
independent of the angle a of the cake shape Andreev billiards. From Fig. Elone can see 
that the length of the longest possible chord is 2R. Therefore, the cut-off of the path 
length distribution P(s) is always 2R and obviously independent of a. 

It is easy to express the gap E gSbP in units of the mean level spacing Sn of the 
isolated normal region. It is well-known that 5n = 2nh 2 /(m^A), where A = ^R 2 a is 
the area of the normal region |22j. Then, from (|27J) one finds 



E gSLp _ ak F R 



(28) 



5 N 16 

for £ c C ii. In a macroscopic sample k F R ~ R/X F ^> 1, where Af is the de Broglie 
wavelength. Thus, the gap in the cake shape Andreev billiards can be a large value on 
the energy scale of the mean level spacing. 
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4. NS interface with potential barrier 

In this section we investigate the dependence of the gap on the potential barrier at 
the NS interface. Here we assume that the effective masses and the Fermi energies 
between the normal and the superconducting regions are the same. Following Mortensen 
et al. [SI] the mismatch of the N and S regions is characterized by the ratio of the 
Fermi wavenumbers and the Fermi velocities of the two regions — fcp ^/^f ^ cind 
r v = v 

rM s) - it 

is easy to show that the ratio of the effective masses and the Fermi 
energies appearing in the secular equation Q can be expressed by the parameters r k 
and r v in the following way ms/m^ = r v /r k and E^ /Ep = r k r v . The energy levels 
of the systems is calculated by solving the secular equation (jHJ) for different Z, and for 
r k = 1 and r v = 1. The result is shown in Fig. |3J It is seen from the figure that the 

0.15 



o.i 

0.05 





12 3 

Z 

Figure 4. Numerically calculated energy gap (in units of A) as a function of Z. 
For Z = the semi-classical result from Q27JI is indicated by the dashed curve. The 
parameters are the same as in Fig. 03 and r k ~\ and r v = 1. 

gap is decreased with increasing potential barrier. This is owing to the fact that the 
probability of the normal reflection is increased and the condition for perfect Andreev 
reflection is supressed. 

5. NS interface in case of mismatch 

We now present our numerical results for the gap when there is a mismatch in the 
effective masses and the Fermi energies of the N and S regions. The numerical solution 
of the secular equation gives the energy levels of the system. The lowest level can 
be associated to the gap. In Fig. EJthe gap is plotted as a function of r k and r v . The 
surface plot shows that the suppression of the gap is very sensitive to the parameter r v 
but practically independent of r k . 
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6. Conclusions 

We calculated the Andreev bound states of cake shape Andreev billiards by solving 
the Bogoliubov-de Gennes equation. Matching the wave functions in the normal and 
superconducting regions we obtained a secular equation for the energy levels of the 
system. In the exact calculation the different effective masses and Fermi energies of the 
two regions, and the potential barrier at the NS interface are taken into account. For 
a perfect interface we derived a semi-classical expression for the step function and the 
DOS of the energy spectrum. These quantities are expressed by the classical return 
probability P(s) of the particles. It was found that the cake shape Andreev billiards 
always possesses a gap. An analytical expression was derived for the value of the gap 
for perfect interface. We also showed that the DOS has singularities at some energies 
given by a simple formula. Since the cut-off of P(s) is the same at which P(s) is 
singular, we found that at the edge of the gap the DOS is singular. Similarly, the gap is 
independent of the angle a of the cake shape Andreev billiard. We think that these facts 
can be exploited experimentally to measure the value of the gap. The Andreev bound 
states determine the tunneling conductance of the NS systems. Thus, the measured 
conductance (proportional to the DOS) should be changed abruptly at the edge of the 
gap. A cake shape Andreev billiard with a = tc can be made by cutting a semi-circle 
region from a semi-infinite bulk superconductor and replacing that by a normal region. 

We also investigated that how the non-perfect NS interface and mismatch in the 
material parameters of the normal and superconducting regions can affect the value 
of the gap. We found that the gap decreases with increasing potential barrier at the 
NS interface. We also calculated the value of the g function of the ratio r k of 

the Fermi wave numbers and the ratio r v of the Fermi velocities of the normal and 
superconducting regions. In Ref. [33j it has been shown that r k determines the ratio 
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of the sine of the incident and outgoing angle of a particle at the NS interface. This 
implies, and it is indeed supported by our numerical results, that the value of the gap 
is very weakly dependent on the parameter r^. However, our numerics show a strong 
dependence of the gap on the parameter r v : increasing the mismatch in the Fermi 
velocities decreases the value of the gap. 
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